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ABSTRACT 

We present a spectroscopic analysis of the black hole binary Nova Muscae 1991 in 
quiescence using data obtained in 2009 with MagE on the Magellan Clay telescope 
and in 2010 with IMACS on the Magellan Baade telescope at the Las Campanas 
Observatory. Emission from the disc is observed in H a, H /3 and Ca II (A8662). A 
prominent hotspot is observed in the Doppler maps of all three emission lines. The 
existence of this spot establishes ongoing mass transfer from the donor star in 2009- 
2010 and, given its absence in the 1993-1995 observations, demonstrates the presence 
of a variable hotspot in the system. We hnd the radial distance to the hotspot from 
the black hole to be consistent with the circularization radius. Our tomograms are 
suggestive of stream-disc overflow in the system. We also detect possible Ca II (A8662) 
absorption from the donor star. 

Key words: accretion, accretion discs - binaries: close - stars: black holes - stars: 
individual: Nova Muscae 1991 - X-rays: binaries 


1 INTRODUCTION 

Nova Muscae 1991 (GU Mus; hereafter Nova Muscae) 
is a low mass X-ray binary discovered in 1991 as a 
bright transient sou rce by instrumen ts on board the Ginqa 
(|Makino et al.|[l9^ 1 and GRAN AT (|Lund fc Brandtlll99lh 
satellites. It is a soft X-ray transient, a class of objects 
which displays changing levels of soft X-rays and infrequent, 
yet extreme X-ray outbursts. Its mass function (/(M) = 
3.02z b0.Q5 Mrr^, Wu et al. in preparation; /(M) = 3.01zb0.15 
M 0 . IOrosz et al.lll996l j is consistent with a black hole as its 
compact object. 

The optical counterpart o f Nova Muscae was discovered 
by Idella Valle. Jarvis Wes3 (|l99ll j who noted a large in¬ 
crease in brightness to V ~ 13.5 in January 1992. The optical 


* This paper includes data gathered with the 6.5 meter Magellan 
Telescopes located at Las Campanas Observatory, Chile, 
t E-mail: cperis@cfa.harvard.edu; c.peris@neu.edu 


spectra taken in 1992 displayed Balmer, He and N lines su¬ 
perimposed on a blue-peaked continuum which established 
it as the optical counterpart to the X-ray nova. 


An important feature of systems such as Nova Muscae is 
the interaction between the in-falling stream and the accre¬ 
tion disc. This interaction typically occurs in a small region 
also known as the “hotspot” and is observed in these systems 
both as a peak in the orbital light curve at phase ~ 0.8 — 0.9 
(where phase 0.0 is defined as the superior conjunction 
of the compact object) and as a bright spot at the same 
phas e in tomograms (e.g., Kr zeminskil Il965l: IW ood et"al] 


19861. Neilsen. St eeghs &: Vrtilek 2 0081 . ICalyelo et ahf 20091 : 


Peris Vrtilekl | 201 ^ . Lubow &: Shul (|l976h were the first 

to suggest that some of the gaseous matter falling from the 
donor star onto the disc via the stream does not stop at the 
hotspot but flows over and under the dis c. Smooth particle 
hydr o dynamics (SPH) simulations (s ee lArmitage Livi3 
1 19961 : iKunze. Soeith k, HessmanI l200l| j have predicted that 
a significant fraction of the stream can be deflected from the 
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edge of the disc and overflow to smaller rad ii. In some cases 
(ICornelisse et alJboOTl : IPearson et al ]|2^ observations of 
neutron star binaries have provided likely confirmation in 
the form of a tell-tale bright feature at phase ~ 0.5 — 0.75 in 
emission line Doppler maps reconstructed via the Doppler 
tomography technique ([Marsh fc HorndllOSSl b We will ex¬ 
plore a possible stream-disc overflow in the accretion disc of 
Nova Muscae. 

Previous optical tomography of Nova Muscae was 
carried out using data obtained with the 3.9 m Anglo- 
Austra lian Telescope and t he 3.6 m New Technology Tele¬ 
scope ([Casares et al.l Il997f ) . ICasares et al.l did not detect a 
hotspot in their tomogram and suggested this may be evi¬ 
dence for low mass transfer rates in Nova Muscae compared 
to other X-ray transients of comparable orbital period. 

In this paper we present optical tomography of Nova 
Muscae in quiescence using data from the 6.5 m Magellan 
Baade and Clay telescopes. In Section 2 we describe our 
observations and data reduction. In Section 3 we analyse 
the spectrum and present the tomographic results while in 
Section 4 we discuss and interpret our results. In Section 5 
we present our conclusions. 


2 OBSERVATIONS 

Nova Muscae was observed with the Magellan Echelette 
(MagE; [Marshall et al.l l2008h spectrograph on the Magel¬ 
lan Clay telescope in 2009. Data were obtained over the 
nights of April 25-26, 33 spectra from the first night and 39 
from the second night. 72 spectra, each containing 15 orders, 
were collected with individual exposure times of 600 s. A slit 
width of 0.85'' was used throughout. The resultant spectral 
resolution achieved is approximately 1 A at a wavelength 
of 5000 A, corresponding to ~60 km s“^ velocity resolution. 
The data reduction steps (bias subtraction, flat-held correc¬ 
tion, wavelength calibration, and spectral extraction) were 
performed using the Carnegie Observatory pipelinqj. The 
sky-subtraction of the 3 highest wavelength orders of the re¬ 
sulting spectrum could not be performed due to limitations 
of the pipeline. 

A set of 20 Nova Muscae observations was obtained us¬ 
ing the IMACS instrument on the Magellan Baade telescope 
in 2010, over the nights of Eebruary 10, 11, and 12. Spec¬ 
tra were collected with exposure times of 1200 s each, using 
a 600 lines mm“^ grating. Atmospheric conditions allowed 
slits of 0.75" (Eeb. 10 and 11) and 1.5" (Eeb. 12). Spectral 
resolutions of 2.6 A and 5.2 A were achieved with these re¬ 
spective slits at H a, corresponding to velocity resolutions of 
~118 km s“^ and ~236 km s“^. Spectra were reduced using 
standard image reduction routines in IRAE version 2.10. 

Erom examination of the images containing the Nova 
Muscae spectra we noted the presence of extended emission 
at H a distributed along the direction of the slit on both 
sides of the spectrum. Diffuse Galactic H a emission is a 
possible explanation; we have tested this conjecture by de¬ 
riving an approximate flux calibration for the observed dif¬ 
fuse emission, and noting that this compares well with the 
H a surface brightness at the position of Nova Muscae in 

^ http://code.obs.carnegiescience.edu/mage-pipeline 


the maps of lEinkbeinerl (|2003l ). We obtained a value for the 
flux of the extended H a emission by comparing the average 
number of counts observed on both sides of the slit to that 
of a standard star (LTT 4816, observed for 300s with the 
same instrument settings as Nova Muscae; we assumed for 
comparison with the diffuse H a that the flux from the point- 
source standard was diffused by seeing over approximately 
1 arcsec^). The flux of diffuse H a around Nova Muscae was 
found to be 5V^1_2_T1 x 10 ~^^ ergs cm“^ s“^ arcsec“^. 
The maps of lEinkbeinerl (|2003l ) show a value of 1.36 x 10“^® 
ergs cm“^ s“^ arcsec“^ for Galactic diffuse H a emission in 
the same region. 

Eor the convenience of future observers, we have up¬ 
dated the coordinates of Nova Muscae using a 60 s image of 
the held obtained using the Bessel V Alter on IMACS. We 
fitted a world coordinate system (WCS) to the image using 
imwcs version 2.8.7 and the 2MASS catalogue was used as 
a WCS reference. Nova Muscae was found to have a right 
ascension of II h 26 m 26.6 s and a declination of -68 d 
40 m 32.87 s (the error of right ascension corresponds to 
0.4 s while the error of the declination is dominated by the 
2MASS uncertaint}Q of 0.07-0.08 arcsec). A finding chart 
with the source marked has been included in the Appendix. 


In Eigure I, we present 10 successfully sky-subtracted orders 
of the mean MagE spectrum, covering a wavelength range of 
~3300-7300 A. A low-order fit to the continuum has been 
subtracted to highlight line features, of which the Balmer 
lines are dominant. Nova Muscae’s Balmer lines display the 
classic double-peaked profile which is a signature of Doppler 
shifted emission originating from the accretion disc. The He 
I lines (5875 and 6678 A) and Ee II multiplet 42 (4924, 5018 
and 5 169 A;|m oorelll97^. noted previously bv lOrosz et al.l 
(|l996l ). are detected. Orosz et al.l detected these features via 
isolation of the disk component of the Nova Muscae spec¬ 
tra by the subtraction of a K5 V template star spectrum. 
Since we do not isolate the disk component we And these fea¬ 
tures at lower S/N. Additionally, we obs erve the He I line at 
4471 A. The NHI/CHI Bowen blend that lOrosz et al.l (|l996l ) 
observed after the 1992 outburst is very weak in our spectra. 

Ca H emission is detected in higher wavelength orders 
of the MagE spectrum which were not successfully sky- 
subtracted. We present one such emission line of interest 
(A = 8662 A) in Eigure 2 with which we produce a tomo¬ 
gram. It should, however, be interpreted with caution due 
to the presence of unsubtracted sky-lines. 

In analogy to Eigure I, Eigure 3 presents the mean 
IMACS spectrum after continuum subtraction. Due to the 
low S/N we could not distinguish any features suitable for 
tomography apart from the strong H a line. Therefore, we 
present only the data from IMACS/CCD2 which covered 
the wavelength range ~6000-6700 A. 


^ http://spider.ipac.caltech.edu/staff/hlm/2mass/overv/overv.html 


3 ANALYSIS 

3.1 Spectrum of Nova Muscae 


© 2015 RAS, MNRAS QQQ.flllTO] 































Tomography of X-ray Nova Muscae 1991 3 


n I r 


3 


50 - 
0 ^ 
-50 - 


1 ^ ^ ^ ^ n 

H3 Hy 

I I He I 

I 




± 


J_I_L_ 


± 


± 


J_I_L_ 


± 


3000 


3500 


4000 


4500 




Wavelength (A) 


Figure 1. The continuum-subtracted mean spectrum of Nova Muscae, from 10 successfully sky-subtracted orders of the MagE data with 
a total integration time of 720 minutes. The spectrum has not been corrected for telluric absorption. Balmer lines dominate the spectral 
profile. 


3.2 Doppler Imaging 

Standard Doppler tomography constructs the distribution 
of line emission from a binary system in a two dimensional 
velocity space image ([Marsh fc HornelflQSSl l . It uses the line 
emission at each orbital phase to calculate the line strength 
as a function of velocity. Figure 4 shows a schematic which 
makes evident the relationship between velocity coordinates 
and spacial coordinates for a disc-and-star system. 

In order to generate the tomograms, each spectrum was 


binned into equal velocity bins and a linear fit performed to 
the continuum around each line. The continuum was then 
subtracted and the residual spectrum provided as input to 
the tomography software, into one of 60 phase interval bins. 
The tomography software modifies a template starting im¬ 
age (a uniform grid of pixels or a two dimensional Gaussian), 
to fit the data by minimizing its reduced in comparison 
to the data. To break the degeneracy inherent in this ap¬ 
proach, the maximum entropy method (MEM) is employed 
to select the smoothest image during computation. 


© 2015 RAS, MNRAS QQQ.fTHTO] 
















4 Peris et al. 
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Figure 3. The continuum-subtracted mean spectrum of Nova Muscae, created using IMACS/CCD2 data taken on Feb 10. The total 
integration time is 160 minutes. 



Wavelength (A) 

Figure 2. The Ca II (A8662) emission line obtained from an order 
of the MagE spectrum that was affected by sky-lines. 


One shortcoming of Doppler tomography is that it as¬ 
sumes emission from all structures is continuous over the 
full binary orbit. Modulation Doppler tomograph^l^ is an 
extension of standard Doppler tomography that allows the 
line flux from a given point to vary as a function of time 
enabling phase-dependent information carried by the line 
profiles to be explored. We produce additional modulation 
Doppler tomograms for each line. 


^ For comprehensive descriptions of Doppler tomography and 
modulati on tomograp hy, r espective l y, we refer the interested 
reader to iMarshI (1 200 if) and ISteeghsI (l2003f) . 


3.3 The Maps 

The orbital and binary characteristics of Nova Muscae have 
been studied elsewhere, with properties summarized in Ta¬ 
ble 1. Using these values as input, in Figures 5-8, we present 
a triplet of tomograms of Nova Muscae, each having been 
computed for a different spectral line: H a, H /3 and Ca II. 
The four panels of each figure depict the following: On the 
upper left, observed line profiles as a function of phase; on 
the upper right, the profiles predicted from the Doppler map; 
on the lower left, the Doppler map; on the lower right, the 
amplitude modulations that result when constant disc emis¬ 
sion is subtracted from the Doppler map. 

Within the panel that contains the Doppler map, we 
over plot the Roche lobe of the secondary and mark two 
tracks representing the stream and disc velocities as in Fig¬ 
ure [H The lower track shows the ballistic trajectory of the 
gas stream and the upper gives the Keplerian disc veloc¬ 
ity along the gas stream. Each open circle marks O.IRli 
and each point repre sents O.QIRl i where Rli is the inner 
Lagrangian distance (|Marshll2005l h The hotspot due to the 
impact of the accretion stream on the disc is usually po¬ 
sitioned on or in-between these lines. The positions of the 
Roche lobes are set using the radial velocity semi-amplitudes 
Ki and i^ 2 . Since Ki has not been directly measured for 
Nova Muscae, we used the mass ratio {q) and K 2 (see Ta¬ 
ble 1) to obtain a value Ki — 27.7 zb 1.2 km s~^ . The main 
purpose for overlaying these trajectories is to interpret the 
position of the hotspot and its relationship to the secondary 
star’s Roche lobe. We therefore tested this placement by 
shifting Ki within its errors and found the displacement to 
be negligible. 

In Figure O we show the tomographic maps for H a, 
the strongest emission line. The disc can be seen clearly in 
the Doppler m ap on the lowe r left p anel. Contrary to the 
observation by ICasares et a D (Il997ll we observe a promi¬ 
nent hotspot on the ballistic trajectory of the gas stream 
at phase ~0.8 where the accretion stream impacts the disc. 
A crescent-shaped emission region from phases ~ 0.2 — 0.6 
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Table 1. System Parameters of Nova Muscae 


Orbital Period, Porb (days) 
Zero ephemeris, To (days) 
Systemic Velocity, 7 
Inclination, i 
Mass Function, f{M) 

Mass Ratio, q 
K2 


0.43260249 ± 0.00000009 
2454946.79462 ± 0.00042 ^ 
14.24 ± 2.08 km s"! 

54 ± 1.5° ° 

3.02 ±0.05 Mq«’* * 
0.068 ±0.003“ 
406.83 ± 2.22 km s"! 


“ Wu et al. (in preparation) Derived using the MagE speetra used in this paper. 
^ Derived from the spectroscopic ephemeris calculated by Wu et al. 

° ^elino^larrisor^^^dcNamam (^OW) 

* Consistent with the values derived by ^ros^^^a^ (^99^ 
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Figure 4. A schematic of the relationship between velocity coordinates (left) and spatial coordinates (right) in a Doppler map. The 
trajectory of the accretion stream is shown as the arced solid line in both panels, and in both the Roche lobe of the secondary star is 
shown in solid gray. The dashed trajectory in the left panel shows Keplerian velocities along the gas stream. In this representation, the 
star orbits counter-clockwise from the perspective of a fixed observer to the right; the number at the position of the secondary star labels 
its phase. The solid-black circle in the right panel marks the outer edge of the disc, at which the orbital velocity in the disc is lowest. 
The Roche lobe of the black hole primary is represented by a dashed line. 


and an additional bright spot at phase ~0.0 are also evident. 
These features are discussed in Section [H 

The H a emission line in the IMACS spectrum (Figure 
6), which yields a similar tomogram, shows a clear hotspot 
as in the MagE H a map (Figure 5). The bright spot at 
phase ~0.0 is also observed. A region of excess emission cor¬ 
responding to the crescent-shaped emission region manifests 
at phases ~0.25-0.6. 


track, shows its velocity to be closer to that of the Keple¬ 
rian disc along the gas stream than the ballistic trajectory. 
This is expected in a situation where the collision of the 
stream with the disc creates a shock which shifts the veloc¬ 
ity of the inflowing stream toward the velocity of the disc. 
Such interme diate positioning has been observed in other 
systems (e.g., Neilsen et al.ll2008l . iMarsh. Robinson Wood 


11994 ICalvdo et allbOOit 


The Doppler map of H /3 on the lower left hand panel 
of Figure 7 also displays the annular emission expected from 
the disc. A bright hotspot is evident from the data at phase 
~0.9. The location of the hotspot in H /3 closer to the upper 


The circular disc structure in the H jS map extends out 
to higher velocities (inwards to smaller radii) than in H a, 
indicative of H /3 production by material at higher temper¬ 
atures and correspondingly closer to the black hole. The 
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•p = 1.267 


= 2.807 



-1000 0 1000-1000 0 1000 


(km/s) (km/s) 

Figure 5. The H ol maps for the MagE data in 2009. The trailed 
spectra is shown in the top-left panel; the Doppler map is on the 
bottom-left; the modulation-amplitude map is on the bottom- 
right; the predicted trail from the Doppler map is on the top- 
right. 


Doppler map shows no other features apart from the hotspot 
and the disc emission. Fluctuation of the hotspot emission 
harmonic to the orbital period is observed in the modulation 
amplitude map. 

Lastly, we present the maps of Ca II (A8662) in Figure 
8 derived from the emission line shown in Figure 2. As men¬ 
tioned in Section 3.1 and evidenced by the affected intensity 
scale, this image should be interpreted with caution due to 
unsubtracted sky lines. However, the presence of clean S- 
curves generated by the bright spots in the predicted trail 
demonstrates that the features in our Doppler map are rela¬ 
tively unaffected by the stationary sky lines. Ca II seems to 
originate from disc material at lower velocities (larger radii) 
when compared to the H a map. This indicates that the 
Ca II emission occurs closer to the outer edge of the disc, 
where one expects lower gas temperatures. The Doppler map 
in Ca II shows the hotspot at phase ~0.85. As in the case 
of H /3, the hotspot in Ca II is offset from the ballistic tra¬ 
jectory of the gas stream. We note the presence of bright 
emission regions near phases ~0.25 and ~0.6 which corre¬ 
spond to the edges of the crescent-shaped emission region 
observed in H a. The hotspot and the other bright regions 
in Ca II lag in phase by ~0.05 when compared to the spots in 
the H a map. We also note the presence of a compact region 
coincident with position of the secondary, which displays 
substantially lower emission relative to the surrounding re¬ 
gion. Fluctuation of the emission from a similar location is 
observed in the modulation amplitude map. 



-1000 0 1000-1000 0 1000 


(km/s) (km/s) 

Figure 6. H a maps created using IMACS data taken in 2010. 
Panels are the same as described in Figure 5. 


^ = 1.000 



-1000 0 1000-1000 0 1000 


(km/s) (km/s) 

Figure 7. The same as Fig. [S] but for H /3, again for the MagE 
2009 spectral data. 

4 DISCUSSION 

4.1 The hotspot of Nova Muscae 

Our tomograms display a bright spot which corresponds to 
the stream-disc interface. The presence of this feature runs 
counter to the hypothesis that Nova Muscae harbours an ab¬ 
normally low rate of mass transfer compared to other black 
hole transients of comparable orbital period, and demon¬ 
strates the presence of a variable hotspot. 

By considering the position of the hotspot relative to 
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X® = 1.098 



-1000 0 1000-1000 0 1000 


(km/s) (km/s) 

Figure 8. The same as Figs. [5] and 7, but for Ca II (8662 A). 
The tomogram was created using an order of the MagE spectrum 
that was affected by sky-lines. 


Figure 9. The H /3 tomogram of Nova Muscae, shown close up. 
The black circle marks EcircSin(i). 



Vcirc — 


GM 


Rc\ 


( 1 ) 


and radius i^circ, given bv I Warn^ (|2003l l for 0.05 <q<l 


the stream’s trajectory, we identify the radial distance to the 
hotspot from the black hole. From the H a Doppler map, the 
hotspot is found to lie squarely along the ballistic trajectory 
of the stream. The open circles along this line mark 10% 
increments in i^Li, which places the position of the hotspot 
near R = 0.5 zb 0.05i7Li- The H [3 Doppler map depicts 
the hotspot closer to the Keplerian trajectory, indicative of 
the disc’s rotation along the course of the stream. This may 
represent post-shock H {3 emission as the stream is swept 
along into the disc. 

Assuming a Keplerian trajectory for the H /3 hotspot, 
we find that the emission originates from R — 0.45zb0.05i^Li, 
in good agreement with the H a map. In the Ca II Doppler 
map, the hotspot appears sandwiched between the ballistic 
and Keplerian curves. The ballistic trajectory gives an up¬ 
per limit for the hotspot at = 0.55 zb O.Obi^Li whereas the 
Keplerian lower limit corresponds to = 0.4zb0.05i^Li. De¬ 
spite the terminology, rather than envisioning a small “spot” 
of emission, one could instead envision the bright interface 
as a stretched patch reflecting a transition between ballis¬ 
tic and Keplerian trajectories, occurring at some distance 
R ~ O.bi^Li. In this case, the consistent distances suggest 
that the emission is relatively confined, but that different 
zones along this interface are responsible for the different line 
features. These zones, which are variously more or less pulled 
into the disc’s motion, would naturally span the gamut from 
ballistic to Keplerian motion, while still originating from ap¬ 
proximately the same location. For the location of such a 
hotspot, we find i^hotspot = 0.5 zb 0.05i^Li. 

In order to more readily place the location of the 
hotspot into context, we overlay the circularization radius 
of the system on the H (3 Doppler map. This radius corre¬ 
sponds to the disc orbit with the same specific angular mo¬ 
mentum as the gas stream. Accordingly, the orbital velocity 
at the circularization radius is given by 


i?circ « 0.0859og"°-^^®, (2) 

where a is the binary separation. For Nova Muscae, Rdvc = 
(8.64 zb 0.22) X 10^° cm (= 0.44 zb O.Oli^Li), and its Keple¬ 
rian velocity Vdrc = 1001 zb 23 km s“^. Kcirc is marked on 
the H (3 Doppler map in Figure 9. Although the tomogram 
indicates that the hotspot emission originates at smaller ve¬ 
locities than Kcirc (i-e. exterior to i^circ), it is noteworthy 
that it lies within 1 a. Considering t hat the maximum dis¬ 
tance to the disc edge is ~ O.Oi^Li ([Frank. King k, Raind 
I 2 OO 2 I : eq. 5.122), Nova Muscae’s hotspot, at ~ 0.5i^Li, is 
consistent with its circularization radius. 

4.2 Possible Stream-Disc Overflow 

SPH simulations performed by lArmitage &: Livi3 (|l996h 
showed that a significant portion of the gas stream can be de¬ 
flected from the outer edge of the disc and flow over the disc. 
The deflected material loses energy and settles at smaller 
radii approaching the circularization r adius. A later analy - 
sis using SPH at higher resolution by iKunze et al.l (|200lh 
revealed that a bow shock at the hotspot launches gas to 
about 20° — 25° above the disc plane. This gas can produce 
absorption dips in the light curve at orbital phases 0.7-0.9 
in systems with inclinations greater than 65°-70°. The sim¬ 
ulations show the deflected mass settling back onto the disc 
by phase ~ 0.5 independent of the mass ratio considered. In 
the right circumstances, such as during a spike in the mass 
transfer rate, or given a smaller sized disc, the simulations 
also show the overflow stream rebounding off the disc at 
phase 0.5 and producing yet another absorption dip in the 
light curve at phase 0.2. 

Our H a (both IMACS in 2009 and MagE in 2010) and 
Ca II maps show a crescent shaped emission feature origi¬ 
nating from the region stretching between phases ~ 0.2—0.6. 
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Figure 10. A rough schematic of the stream-disc overflow in 
Nova Muscae. The numbers signify the phase of the system rel¬ 
ative to the observer when the donor star reaches the marked 
position. The yellow spots indicate the hotspot where the accre¬ 
tion stream first hits the disc and the point of fallback of the 
deflected flow. The rise in elevation of the deflected flow above 
the disc is indicated by vertical lines. 


However, further investigation is hindered by the presence 
of u nsubtracted sky lines in our data. 

ICasares et al.l (|l997f ) detected emission in H a from the 
region of the donor star, which they attributed to possible 
chromospheric activity. We detect a similar bright spot, but 
find that it is offset from the secondary at AV ~ 200 km s“^ 
larger Vy owing to the updated dynamical model. H a emis¬ 
sion, identified as originating from the secondary, has been 
observed in other b lack hole and neutron star systems (e.g. , 
Marsh et al. 1994: Shahbaz et al.lll999l: lTor res_eHahl 200d : 
Gonzalez Hernandez Casaresll201^ : IPeris Vrtilekll2012l i , 


some compatible with the Roche lobe of the secondary and 
some slightly offset. Although it has been commonly at¬ 
tributed to chromospheric activity or X-ray heating of the 
secondary, an unequivocal interpretation has not yet been 
established. Our observation of similar H a emission with 
a large offset is an intriguing result and we do not have an 
explanation for this as of yet. 


The associated region is consistent with the point at which 
the deflected m ass stream falls back onto the disc as in 
the iKunze et al.l picture. Together with the presence of the 
hotspot, this is suggestive of stream-disc overflow. The shift 
in velocity of the heated stream as evidenced by the hotspot 
drift from F igures 7-8 is consistent with the results of the 
iKunze et al.l simulations as well. 

Within this overflow paradigm, one expects a diffuse 
cloud of absorbing gas to persist above the disc, reaching 
zenith at phase ^0.7. Altho ugh the inclination of Nova Mus¬ 
cae, 54° (|Gelino et al.ll200lh . is likely too low to manifest cor¬ 
responding absorption dips in its light curve, the absorbing 
matter is expected to be evident in H a emission. Our H a 
Doppler map also displays a gap in which the line emission is 
reduced around phase ~0.7 and ~ 0.15, th e phases expe cted 
for absorption from mass overflow in the iKunze et al.l pic¬ 
ture. ft is interesting to note that t h e H a equivalent width 
variations shown by iGasares et al.l (|l997l i in their figure 3 
produce dips at phases ~0.7 and Their and our ob¬ 

servations are in good accord with the iKunze et al.l simula¬ 
tions. 

In Figure 10 we present a rough schematic of the over¬ 
flow process. We note that because velocity coordinates do 
not uniquely map into spatial coordinates, this figure is only 
approximate, based on one plausible realization of our tomo¬ 
grams. 

4.3 Detection of the Secondary Star 

While the Ga 11 emission from the disc is clearly present, 
our Ga 11 Doppler map also displays a compact region of 
lower emission coincident with the position of the secondary. 
We also observe fluctuation of the flux originating from a 
similar position in the modulation amplitude map. This is 
consistent with Ga 11 absorption from the secondary and 
is strongly reminiscent of the similar Ga 11 emission an d 
absorption seen in XTE J1118+480 bv iGalvelo et al.l (|2QQ9li . 

^ After accounting for a phase shift of 0.06 to bring their data 
up to date. 


4.4 Comparison between Nova Muscae and 
A0620-00 

In this se ction we compare our tomograms wit h those of 
A0620-00 (|Neilsen et ahllioosl : iMarsh et al.lll994l h a system 
which is widely compared to Nova Muscae in the literature 
given the similarity of the system characteristics. The orbital 
period of Nova Muscae (Ppr b — 10.38h) slightly exceeds tha t 
of A0620-00 (Porb = 7.75h: [McGlintock Remillardlll986l i. 
When in outburst, A0620-00 has been suggested to harbour 
a slightly thicker and hotter disc compared to Nova Muscae 
(|Esiii et al 

Since both systems display prominent hotspots in 
Balmer emission, it is straightforward to compare them us¬ 
ing the Doppler maps as a guide. To keep the comparison 
robust, we focus our comparison on examining the positions 
of the hotspots for both systems. We e xamine tomograms 
in H a, figure 10 in iNeilsen et al.l (|2008h and Eigure 5 & 6 
here. As described in § 4.1, we find a hotspot radius for Nova 
Muscae of R — 0.5 ± 0.05Rli; for A0620-00, the H a map 
reveals a larger but consistent radius, R = 0.6 ± 0.05Rli. 
Assuming that the H [3 emission originates from circular¬ 
ized gas in Keplerian orbit. Nova Muscae’s hotspot corre¬ 
sponds to a radius of R = 0.45 ± 0.05Rli and A0620-00’s 
hotspot to R = 0.5 ± 0.05Rli. These results are suggestive 
of a marginally smaller distance to Nova Muscae’s hotspot, 
relative to the binary separation. 

SPH simulations conducted by lEoulkes et al.l (|2004h 
showed that eccentric accretion discs manifest several char¬ 
acteristics in a Doppler tomogram. Notably, as the asymmet¬ 
ric disc rotates, the stream-disc impact point travels along a 
section of the edge of the disc creating an extended hotspot. 
Two distinct non-sinusoidal S-curves are then created by the 
disc gas and stream gas emission at the hotspot. They also 
predict crescent-shaped emission from the disc in the tomo¬ 
gram. An eccentric accretion disc should also show an offset 
of its centre of symmetry from the central position of the 
black hole. 

Although the H a tomogram in Nova Muscae displays 
a crescent shaped region of emission, we see no other com¬ 
pellin g evidence for disc eccentricity. However. iNeilsen et al.l 
(|2008l i finds compelling tomographic evidence that the disc 
in A0620-00 is eccentric. We suggest that unlike in the case 
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of A0620-00 llNeilsen et, alJlioO^ the disc of Nova Muscae 
might not have grown large enough for the growth of eccen¬ 
tric modes. 


5 CONCLUSIONS 

Our high resolution tomography of Nova Muscae shows a 
hotspot in both 2009 and 2010, confirming ongoing mass 
transfer via Roche-lobe overflow in the binary system during 
this epoch. Given its earlier absence, this raises the interest¬ 
ing possibility of a variable hotspot in Nova Muscae. 

The Ca II maps display evidence for absorption from 
the secondary as well as emission from the disc. A s imilar 
observation for XTE J1118+480 (jCalvelo et al.|[2009l l sug¬ 
gests that Ca II (A8662) is of particularly high utility in 
cons training binary parameters in X-ray binary systems (see 
also Ivan Spaandonk et al.l l2010l ) . Higher resolution spectra 
of the Ca II will help to better map the secondary in Nova 
Muscae. 

Our tomograms support stream-disc overflow in the sys¬ 
tem. The bright spot in the lower half of the H a map (also 
found in approximately the same location in the Ca II map) 
at phase ~0.6 is interpreted as the point at which the de¬ 
flected mass falls back to the disc. The overflow should re¬ 
sult in an elevated cloud of absorbing gas above the disc, 
reaching its maximum altitudes at phase ~0.7. Therefore, 
the presence of an absorption dip in the light curve at phase 
~0.7 would provide a “smoking gun” for this interpretation. 
However, due to the moderate inclination of the system, the 
lack of a dip would not preclude this scenario. Such over¬ 
flow might also ex plain the broader blue wing in the H a 
profile observed bv lOrosz et al.l (|l992l i which was previously 
interpreted as an outflow of mass from the disc. 

Bright spots in the phase region ~0.5 - 0.75 have often 
been observe d in other systems in Balmer and He I lines 
Jm ars i] l200lh . We assert that these may be due to the ubiq¬ 
uitous presence of stream-disc overflow causing heating of 
the surface of the disc in the re-impact region. 


ACKNOWLEDGEMENTS 

We gratefuly acknowledge the use of the MOLLY software 
package written by Tom Marsh. We would like to thank 
Jorge Casares for sharing his data on Nova Muscae with 
us. We would also like to thank the referee for helpful com¬ 
ments. JFS has been supported by NASA Hubble Fellow¬ 
ship grant HST-HF-5I3I5.0I. LCH acknowledges support 
from the Kavli Foundation, Peking University, and the Chi¬ 
nese Academy of Science through grant No. XDB09030I02 
(Emergence of Cosmological Structures) from the Strategic 
Priority Research Program. 


REFERENCES 

Armitage P. J., Livio M., 1996, ApJ, 470, 1024 
Calvelo D. E., Vrtilek S. D., Steeghs D., Torres M. A. P., 
Neilsen J., Filippenko A. V., Gonzalez Hernandez J. L, 
2009, MNRAS, 399, 539 

Casares J., Martin E. L., Charles P. A., Molaro P., Rebolo 
R., 1997, New Astron., I, 299 


Cornelisse R., Steeghs D., Casares J., Charles P. A., Barnes 
A. D., Hynes R. L, O’Brien K., 2007, MNRAS, 380, 1219 
della Valle M., Jarvis B. J., West R. M., 1991, Nature, 353, 
50 

Esin A. A., Kuulkers E., McClintock J. E., Narayan R., 
2000, ApJ, 532, 1069 
Finkbeiner D. R, 2003, ApJS, 146, 407 
Foulkes S. B., Haswell C. A., Murray J. R., Rolfe D. J., 
2004, MNRAS, 349, 1179 

Frank J., King A., Raine D. J., 2002, Accretion Power in 
Astrophysics: Third Edition 

Gelino D. M., Harrison T. E., McNamara B. J., 2001, AJ, 
122, 971 

Gonzalez Hernandez J. L, Casares J., 2010, AVA, 516, A58 
Krzeminski W., 1965, ApJ, 142, 1051 
Kunze S., Speith R., Hessman F. V., 2001, MNRAS, 322, 
499 

Lubow S. H., Shu F. H., 1976, ApJ, 207, L53 
Lund N., Brandt S., 1991, lAU Circ., 5161, I 
Makino F., McNaught R. H., Jones A., West R. M., 1991, 
lAU Circ., 5161, I 

Marsh T. R., 2001, in Boffin H. M. J., Steeghs D., Cuypers 
J., eds, LNP Vol. 573, Astrotomography: Indirect Imaging 
Methods in Observational Astronomy, Springer-Verlag, 
Berlin, p. I 

Marsh T. R., 2005, ApVSS, 296, 403 
Marsh T. R., Horne K., 1988, MNRAS, 235, 269 
Marsh T. R., Robinson E. L., Wood J. H., 1994, MNRAS, 
266, 137 

Marshall J. L., Buries S., Thompson LB., Shectman S. A., 
Bigelow B. C., Burley G., Birk C., Estrada J., Jones P., 
Smith M., Kowal V., Castillo J., Storts R., Ortiz G., 2008, 
in SPIE Conference Series Vol. 7014, 701454 
McClintock J. E., Remillard R. A., 1986, ApJ, 308, 110 
Moore C. E., 1972, A Multiplet Table of Astrophysical In¬ 
terest. National Bureau of Standards, Washington, p. 66 
Neilsen J., Steeghs D., Vrtilek S. D., 2008, MNRAS, 384, 
849 

Orosz J. A., Bailyn C. D., McClintock J. E., Remillard 
R. A., 1996, ApJ, 468, 380 

Orosz J. A., Bailyn C. D., Remillard R. A., McClintock 
J. E., 1992, in Bulletin of the American Astronomical So¬ 
ciety Vol. 24, p. 1236 

Pearson K. J., Hynes R. L, Steeghs D., Jonker P. G., 
Haswell C. A., King A. R., O’Brien K., Nelemans G., 
Mendez M., 2006, ApJ, 648, 1169 
Peris C. S., Vrtilek S. D., 2012, MNRAS, 427, 1043 
Shahbaz T., van der Hooft F., Casares J., Charles P. A., 
van Paradijs J., 1999, MNRAS, 306, 89 
Steeghs D., 2003, MNRAS, 344, 448 

Torres M. A. P., Casares J., Martmez-Pais 1. G., Charles 
P. A., 2002, MNRAS, 334, 233 
van Spaandonk L., Steeghs D., Marsh T. R., Torres 
M. A. R, 2010, MNRAS, 401, 1857 
Warner B., 2003, Cataclysmic Variable Stars 
Wood J., Horne K., Berriman G., Wade R., O’Donoghue 
D., Warner B., 1986, MNRAS, 219, 629 


© 2015 RAS, MNRAS QQQ.mfTOl 













10 Peris et al. 



”o.O 35.0 11:26:30.0 25.0 20.0 15.0 10.0 


Right ascension 


Figure 11. The field of Nova Muscae 1991 with the system marked by two lines. 


APPENDIX 

Figure 11 shows a finding chart of Nova Muscae created 
using a 60 s image obtained using the Bessel V filter on 
IMACS (~0.85 arcsec seeing). 
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